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Quantum plasmonics relies on a new paradigm for light-matter interaction. It benefits from strong 
confinement of surface plasmon polaritons (SPP) that ensures efficient coupling at a deep subwavelength 
scale, instead of working with a long lifetime cavity polariton that increases the duration of interaction. The 
large bandwidth and the strong confinement of one dimensional SPP enable controlled manipulation of a 
nearby quantum emitter. This paves the way to ultrafast nanooptical devices. However, the large SPP 
bandwidth originates from strong losses so that a clear understanding of the coupling process is needed. In 
this report, we investigate in details the coupling between a single emitter and a plasmonic nanowire, but 
also SPP mediated coupling between two emitters. We notably clarify the role of losses in the Purcell factor, 
unavoidable to achieve nanoscale confinement down to 10 ~ (Mri) 3 . Both the retarded and band-edge 
quasi- static regimes are discussed. 
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Light matter interaction at the single atom/photon level is extremely weak and strategies have to be developed 
to enhance this coupling. In a simple picture, the coupling strength depends on the quality factor Q and the 
effective volume V of the involved mode by the ratio Q/V. High Q/V ratio, hence efficient light-matter 
interaction, is achieved in cavity-based systems, where large modal volumes are easily compensated by extremely 
narrow resonances 1 ' 2 . It has been also proposed to work with diffraction limited volumes by coupling a single- 
atom to an elongated fiber 3 . Another downscaling step is made with plasmonics 4 . Surface plasmon polariton 
(SPP) results from coupling electromagnetic wave to a surface density of charges. They are therefore naturally 
confined near a metal surface and are not diffraction limited, offering thus a new tool for interfacing light and 
matter at the nanometer scale 5 " 7 . Since the seminal work of Chang et al 8 , coupling between a dipolar emitter and a 
metal nanowire received strong interest. Metal nanowires define ID plasmonic waveguides with a great potential 
for integrated optical routing 9 . They can be chemically synthesized with high crystallinity hence supporting SPP 
with reduced losses 10 " 12 . Controlled positionning of quantum emitter(s) near ID plasmonic guide is now achiev- 
able 13 " 15 , opening the way to realize original nano-optical devices. For instance, SPP propagation can be controlled 
by a single emitter, leading to the concept of single photon transistor 16 . Reciprocally, two distant emitters can be 
interfaced via surface plasmons, with applications such as SPP mediated resonant energy transfer 17 , remote qubits 
entanglement 1819 or nano-optical logical gates 20 . 

This report is devoted to a careful analysis of the coupling mechanisms between one or two individual emitters 
and a metal nanowire. We quantify all the coupling channels (radiative and non radiative coupling, surface 
plasmon coupled emission, electron scattering). To date, radiative and non radiative rates were estimated within 
the quasistatic approximation 8 . Moreover, the expression of the decay rate into plasmon modes were derived in 
the lossless case only 8,21 . In this context, we are specifically interested in the role of losses in the emitter-SPP 
coupling mechanism. Therefore, we investigate two distinct regimes. A first regime taking into account retarda- 
tion explaining SPP propagation and routing, and a second regime focused on quasi-static SPP responsible for an 
efficient interfacing with dipolar emitters but at the price of strong losses. To this aim, we consider a gold 
nanowire since both retardation regime and quasi- static pictures are present in the optical domain. Figure 1 
represents the studied configuration. A dipolar emitter excites a surface plasmon supported by a gold nanowire. 
The report is organized as follows. We first carefully describe the modes properties of a ID plasmonic nanowire in 
presence of losses. This allows for a rigorous description of the coupled emitter-nanowire system in the retarded 
regime. We identify all the relaxation channels of the emitter close to ID plasmonic nanowire and compare our 
data to previous works. We notably achieve a closed form expression for surface plasmon coupled emission, in 
presence of losses with a clear link to the Purcell factor. This closed form expression permits to unambiguously 
define the mode volume of /ossy ID SPP. Then, we discuss SPP mediated dipole-dipole coupling between two 
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Figure 1 | Surface plasmon coupled emission of a dipolar quantum 
emitter close to a metal nanowire of radius R. We represent the intensity at 
a given instant of the guided SPP. p refers to the emitter dipole moment 
and d is the distance to the nanowire surface. €\ and e 2 refer to the dielectric 
constant of the surroundings and the metal nanowire, respectively. 

emitters. In the last section, we fix the operating at k = 525 nm, a 
wavelength corresponding to the quasi- static regime at the band edge 
of the dispersion relation. 

Results 

In order to discuss the relaxation channels available for coupling the 
emitter into the plasmonic waveguide, we first introduce the concept 
of local density of modes in link with the density of modes supported 
by the metal nanowire. 

Metal nanowire: modes density, retarded and quasi-static regimes. 

Figure 2 represents the dispersion relation of the TM 01 mode 
supported by a gold nanowire. Strictly speaking, the stationary 
mode is defined in absence of dissipation (infinite lifetime), as 
shown in Fig. 2a. This surface plasmon polariton results from 
coupling a photon and a surface plasmon so that its dispersion rela- 
tion closely follows the light line at low frequencies (propagating SPP, 
retarded regime) and bends towards the quasi-static SPP resonance at 
larger frequencies. Indeed, in the quasi-static approximation, the SPP 
resonance condition obeys e 2 (co) +£\ = 0, satisfied for co = 3.5 
10 15 rad.s~ l (k = 525 nm), in the present case. The SPP presents a 
zero group velocity in this regime (in absence of loss). The density of 
guided modes (TM 01 ), proportional to dkjdco, is represented in 
Fig. 2b. It reveals a strong accumulation of modes near co = 3.5 
10 15 rad.s~ l that defines a common domain with the high-order 
localized quasi-static modes sustained by a metal nanoparticle. In 
the following, the quasi-static regime refers to the SPP behaviour 
near k = 525 nm. 



In realistic conditions, the SPP mode must be defined in presence 
of absorption as shown in Fig. 2c. This defines a quasi mode with a 
finite lifetime. We observe a similar dispersion curve except near the 
quasi-static asymptote where a typical back-bending occurs 23 . This 
area corresponds to a mode lifetime shorter than its oscillation per- 
iod, or equivalently, a propagation length shorter than the spatial 
oscillation (see Fig. 2d). The concept of modes is then unavailable. 
However, this part of the dispersion curve is essential to describe non 
radiative coupling of a dipolar emitter to a metal nanostructure. We 
therefore abusively associate the term SPP mode to the whole dis- 
persion curve. The retarded and quasi- static regimes are of great 
interest to guide (propagating SPP) or confine (quasi-static SPP) 
light at the nanoscale. The mode propagation length is finally repre- 
sented in Fig. 2d). As expected, the quasi-static SPP is strongly con- 
fined on both the transverse [high k z , hence short penetration depth 

in air dspp~l j (k 2 z — e 2 kiy^ 2 ] and longitudinal (low propagation 

length L SPP ) directions. An extremely small modal volume can thus 
be defined so that efficient coupling to a nearby quantum emitter is 

expected 8,24 . 

Local density of modes, Purcell factor and relaxation channels 
near a plasmonic wire. The local density of modes (2D-LDOS) 
near a plasmonic or optical waveguide reads 25 

Ap 2D (r h k z ) = ~ ^ImTr£(rn)AG 2D (r||,r||,^), (1) 

where AG 2D refers to the Green's dyad part associated with the 
waveguide structure only (i-e excluding its surroundings), c(r\\) is 
the dielectric constant at the location r\\ =R + d and k z is the 
wavevector along the nanowire axis. This definition allows to 
evaluate the density of guided modes (that cancels for k z = 0) but 
also includes all the scattering and non radiative channels that are 
excited by a nearby dipolar emitter. This is therefore extremely useful 
to quantify and describe the dipolar emitter relaxation. Figure 3a 
shows the calculated 2D-LDOS near a gold nanowire as a function 
of its radius at the wavelength k = 800 nm. We observe two 
resonances corresponding to the TM 01 and HE U modes, 
respectively. A strong increase of the (TM 01 ) 2D-LDOS appears for 
small radii that corresponds to a quasi -static -like behavior (different 
from the quasi- static regime that occurs at a specific wavelength and 
investigated later). The mode characteristics (wave vector and 
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Figure 2 | (a) Dispersion relation co = j{k z ) of a gold nanowire of radius 100 nm in absence of losses. k z is the propagation constant of the guided 
mode at an angular frequency co. The horizontal asymptote at the angular frequency co = 3.7 X 10 15 rad.s~ l corresponds to e 2 (co) +£i =0 (quasi-static 
regime), (b) DOS of the TM 01 mode deduced from dkjdco (in absence of loss), (c) Dispersion relation taking into account losses, (d) Comparison 
of the TM 01 mode propagation length (L SPP ) and spatial period (k S pp). The mode propagation length is shorter than its oscillation period for 
co > 3.10 15 rad.s~ l (grey area). The gold dielectric constant is taken from Ref. 22. The dielectric constant of the surrounding medium is €\ = 2.25. 
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Figure 3 | (a) 2D-LDOS calculated at distance R/2 from the surface of a gold nanowire as a function of its radius R. The superimposed white lines 
represents the dispersion relations of the waveguided modes (taking into account the losses, (b) 2D-LDOS as a function of the normalized propagation 
constant k z lk Q for different metal absorption Imfa). The nanowire radius is fixed at R = 10 nm. The 2D-LDOS follows a Lorentzian profile peaked at the 
mode effective index n S pp. The FWHM An S pp is inversely proportional to mode propagation length L SPP . (c) Evolution of the ratio Ap 2D /L SPP as a 
function of Imie-i). 



propagation length) are easily extracted from this formalism 25 . 
Indeed, for a given radius, the 2D-LDOS presents a Lorentzian 
profile near SPP resonance. It is peaked at the mode propagation 
constant and with a full- width- at-half-maximum (FWHM) inversely 
proportional to the propagation length as shown in Fig. 3b. As 
expected, the propagation length strongly decreases by increasing 
the metal absorption Im(e 2 ) and the 2D-LDOS profile flattens. 
However, the mode effective index h S pp does not depend on the 
losses. The number of supported modes is given by the area below 
the 2D-LDOS curve. For a Lorentzian shape, it is proportional 
to the ratio Ap 2D /L SPP . This ratio is therefore independent on 
absorption over a very wide range as observed in Fig. 3c. For large 
losses, the 2D-LDOS does not follows a Lorentzian profile anymore 
and the concept of mode is then unavailable. In the non- absorbing 
limit, the 2D-LDOS profile associated to the SPP resonance is a Dirac 
distribution. 

Following Fermi's golden rule, we can now quantify the modifica- 
tion of the fluorescence rate for a point-like dipolar emitter coupled 
to the nanowire. Indeed, its spontaneous emission is proportional to 
the 3D-LDOS at the emitter position. The 3D-LDOS is simply related 
to the /c z -Fourier transform of the 2D-LDOS so that the total decay 
rate is achieved by integrating the 2D-LDOS on the whole wave- 
vector spectrum (see methods). Physically, the region \k z \ < riik 0 
corresponds to propagative waves into the medium 1, so-called radi- 
ative waves whereas \k z \ > riik 0 corresponds to surface waves prop- 
agating along the nanowire (including SPP and so-called lossy 
surfaces waves) 27 . The 2D-LDOS beyond \k z \ > riik 0 contributes to 
the non-radiative rate only but a small amount of radiative waves 
(\k z \ < riik 0 ) could also be absorbed into the metal nanowire. The 
radiated power is given by the Poynting vector flux. Its computation 
thanks to the Green's dyad technic involves an integration over the 
full nanowire volume (or at least on several propagation lengths) 28 , 
that is a challenging computing task for the infinite wire considered 
here. Therefore, the radiative rate is evaluated by integrating the 2D- 
LDOS on the radiative spectrum (\k z \ < nik 0 ) and cancelling the 
metal losses [lm(e 2 ) = 0] and is represented in Fig. 4a. The non- 
radiative rate is the difference between total and radiative rates and 



is represented in Fig. 4b. We also represent the results from a quasi- 
static approximation for comparison 29 . 

Let us first discuss the radiative rate behavior (Fig. 4a). T rad is only 
weakly modified in presence of the metal nanowire. The quasi-static 
approximation reproduces very well the distance dependence, except 
for a dipole parallel to the wire axis (Fig. 4a,iii) but the modification 
remains low. Indeed, quasi- static approximation considers the dipole 
induced in the nanowire that cancels along the longitudinal axis 29 . 
The dependence of radiative decay rate with distance closely resem- 
bles to the radiative decay rate near a spherical metal particle 30 for 
both radial (Fig. 4a,i) and orthoradial (Fig. 4a,ii) dipoles since the 
decay rate is governed by the dipole induced in the metal nanopar- 
ticle for both geometries. 

The non radiative rate is presented in Fig. 4b. Two main channels 
contribute to the non radiative coupling 26 ' 27 : excitation of the guided 
SPP with a rate T SPP and lossy surface waves at large k z (see also 
Fig. 3a at small radii where high a 2D-LDOS appears for k z = k SPP 
and above). Although included in the non radiative contribution, the 
surface plasmon coupled emission exists even for lossless metal and 
is related to the excitation of a guided mode that can be detected at 
the nanowire output. This point is discussed in details later. In the 
case of a realistic absorbing metal, the SPP part contributes to the non 
radiative channel. SPP absorption originates from the finite conduc- 
tivity of the metal. It dissipates the energy of the current density 
associated to the SPP propagation. Therefore, it is referred as Joule 
(or ohmic) losses. The second dissipation process appearing at large 
momenta is electron scattering at the lattice impurities. Finally, elec- 
tron-hole excitation in the metal occurs at very short distances 27 . 
Note that lossy surface waves describe the dissipated power from 
both electron scattering and electron-hole excitation. The modelisa- 
tion of electron-hole excitation would necessitate to include non- 
local effects in the permittivity. This process is negligible for separa- 
tion distances above d~5 nm considered here 27 . 

At very short distances, electron scattering loss dominates (r SPP 
<C T NRad ). Above 10 nm, the non radiative channel is dominated by 
coupling into the propagating SPP (Fig. 4b). We also observe a small 
contribution of the HE U mode to T Nrad above 80 nm in case of 



SCIENTIFIC REPORTS | 3 : 2734 | DOI: 1 0.1 038/srep02734 



3 



3 




i— ^ 

3 



0 20 40 60 80 100 

d (nm) 
b) i) 




1.05 



■h 0.95 



L 0 20 40 60 80 100 

J (nm) 



10 z 
10 1 
10° 



u io- 1 



10- 



\ '' lin< ; ar 2oo 








* scale 








\ 100 


1 






\ 


0 5 10 15 


\ \ 




\ 




\ 




s 




S *v 



















i 



20 40 60 80 

d (nm) 



100 




20 40 60 80 

d (nm) 



100 



0.85 




0 20 40 60 80 100 

d (nm) 




0 20 40 60 80 

d (nm) 



Figure 4 | Normalized radiative (a) and non-radiative (b) rates of a dipolar emitter placed at a distance d to a gold nanowire for three emitter 
orientations. The nanowire radius is R = 10 nm and the emission wavelength is k = 800 nm. The surface plasmon coupled emission contribution to the 
non radiative rate is also indicated. For orthoradially oriented emitter r SPP = 0 [Fig. b,ii)], since guided mode does not possess orthoradial electric 
component. QS refers to the quasi-static approximation. 



radially oriented dipole [see the small difference between T Nra d and 
T S pp in Fig. 4b(i) and Fig. S4 in supplementary information]. For 
radially (Fig. 4b,i) or longitudinallly (Fig. 4b,iii) oriented emitter, the 
quasi-static approximation strongly underestimates the non radi- 
ative rate since it does not include SPP contribution. In the case of 
an orthoradial emitter, the SPP contribution is null and the quasi- 
static approximation is slightly better (Fig. 4b,ii). 

The surface plasmon coupled emission finally expresses as 26 

T SPP = 3nk Ap 2 u D (r h k S pp) 
HxTq 4n\kspp L SPP 

where both the forward and backward propagation modes are 
included. This expression extends the definition of the Purcell factor 
to a guided mode picture making a useful link with cavity quantum 
electrodynamics. As expected, the coupling rate depends on the 
overlap with the mode profile (represented by the density of guided 
mode Ap 2D ) so that it decreases with the emitter-nanowire distance 
(penetration depth in air at k = 800 nm: d SPP ~ 26 nm). The l/L SPP 
dependence is more surprising at first glance. Close inspection of the 
expression clarifies this point. This Purcell factor (Eq. 2) is constant 
on a large absorption range, since the ratio Ap 2D /L SPP determines the 
number of modes and does not depend on the losses (see Fig. 3b,c). 
This implies that an emitter couples to the waveguide independently 
on the following dissipations: Joule losses or propagation up to the 
waveguide output 26 . Obviously, lower absorption is preferred to 



achieve a sufficient signal level at the waveguide output. This original 
result can also be understood from the expression of the Purcell 
factor 

f = rspp = 3 f x y ° as 

where Q and V e ff refer to the mode quality factor and effective 
volume, respectively. Here, the mode quality factor is Q = k SPP l 
Ak SPP = n SPP k 0 L SPP (typically Q^^for.R^lOnm and Q = 62 
for R = 100 nm, see supplementary information). The mode volume 
is estimated as follows. The mode confinement in the transverse 
plane is characterized by an effective surface A ej ^that we explicit later 
on. Due to the finite propagation length L SPP in the longitudinal 
direction, the mode volume varies as V e ff ~ A e jj X L SPP . Finally, 
the Purcell factor F p oc QIV e ff ^ A e jj does not depend on L SPP . 
Identifying expressions (2) and (3), we can deduce the mode volume 
of the guided SPP. For this, we assume a randomly oriented emitter 
located at the nanowire surface (see also Ref. 24). Since the Purcell 
factor does not depend on Joule losses, the mode area is computed 
using the lossless coupling rate description (see supplementary 
information). We represent in Fig. 5a the guided SPP mode volume 
for nanowire radius varying from 10 nm to 100 nm. Its is signifi- 
cantly subwavelength demonstrating the ability of surface plasmon 
to confine light at the nanoscale. It is however worthwhile to note 
that the mode confinement along the longitudinal direction is assoc- 
iated to the finite propagation length L SPP . Since it is compensated by 
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the mode quality factor Q in the Purcell factor, the only parameter of 
interest is the mode effective area represented in the inset of Fig. 5a. 
As a partial conclusion, we stress that the plasmonic Purcell factor 
has to be manipulated with care and cannot be used as a figure of 
merit of the coupled system without an additionnal analysis of losses. 
Figure 5b presents the mode confinement at k = 525 nm for com- 
parison and will be discussed in the last paragraph. 

Dipole-SPP coupling efficiency and SPP mediated dipole-dipole 
coupling. Thus far we investigated in details the various relaxations 
channels for a dipolar emitter coupled to a metal nanowire. Finally, a 
key parameter is the ^-factor that quantifies the coupling efficiency 
into the guided SPP as compared to all available channels. /? = r SPP / 
(^Nrad + ^rad) is represented in Fig. 6a. At very short distances 
electron scattering loss dominates. At long distances the plasmonic 
decay rate becomes negligible compared to the radiative rate. But for 
emitter-nanowire distances in the 10-30 nm range, a high coupling 
efficiency above 85% is achieved for both longitudinal and radial 
orientations. As discussed previously, surface plasmon coupled 
emission T SP p does not depend on Joule losses. Nevertheless, low 
loss systems such as crystalline nanowire are evidently preferred. We 
recently demonstrated that the plasmon coupling rate to a pentat- 
winned crystalline silver nanowire is optimized near the nanowire 
edges 12 ' 26 . This means that the plasmon coupling rate r SPP is similar, 
or even better, to the case studied here. Additionally, crystalline 
nanowire presents lower losses rate 10 T NR leading to higher 
factor and signal level at the plasmonic waveguide output. 

Since excellent surface plasmon coupled emission is achieved, we 
now consider two coupled emitters placed near the nanowire at 
distance Az from each other. The coupling rate T 12 between two 
radial emitters is represented in Fig. 6b as a function of the emitter's 
separation distance. If a donor emitter in its excited state is weakly 
coupled to an acceptor emitter in its ground state, T 12 refers to the 
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coupling rate between the two emitters. This weak coupling leads to 
an irreversible resonant energy transfer from the donor to the 
acceptor 3134 . Differently, if the two emitters are strongly coupled, 
the system cannot be described from isolated dipolar emitters in their 
ground (g) or excited (e) state. One rather have to introduce the 
coupled basis \ee>, \gg> 9 | + > = (\eg> + \ge>)/y/(2j and 
I — > = (\ e g> ~ \g e> )/ \f(?) f° r identical emitters. The two 
entangled states | + > and | — > are of particular interest since they 
do not possess classical counterparts. T 12 is the so-called cooperative 
decay rate of the coupled system and r n ± T 12 refer to the decay 
rates of these coupled states. Positive and negative cooperative rates 
r 12 correspond to superradiance and subradiance of the two coupled 
emitters, respectively (note that the physical decay rates r n ± T 12 
are always positives). We represent in Fig. 6b the cooperative decay 
rate. We also evaluated the radiative and SPP contributions to the 
collective emission. The highest collective decay rate (Fig. 6b) corre- 
sponds to the collective excitation of guided SPP but also lossy waves. 
The radiative contribution to the collective relaxation is negligible. 
Therefore the useful part of the collective decay rate is the collective 
surface plasmon coupled emission, achieved for the highest ^-factor, 
that is for d = 22 nm in the present case (Fig. 6a,c). For emitters close 
to the nanowire the coupling rate T 12 presents an oscillatory behavior 
with a spatial period corresponding to the SPP wavelength (Fig. 6b,c). 
For larger separation distances, the cooperative decay rate presents a 
more complex behaviour due to a beating between free- space direct 
emission and the surface plasmon coupled emission (see Fig. 6d and 
S8). 

Quasi-static regime. In this last section, we investigate the surface 
plasmon coupled emission at the band edge of the dispersion 
relation, corresponding to the quasi-static behavior of the SPP (see 



Fig. 2). Indeed, the mode group velocity strongly slows down in this 
spectral range, and we therefore expect important light-matter 
interaction enhancement in this regime 19,35 " 37 . 

Radiative and non radiative rates are calculated in Fig. 7 for an 
emission wavelength X = 525 nm. The radiative rate is significantly 
enhanced near the surface plasmon resonance (Fig. 7a). As expected, 
the quasi- static approximation better describes the non radiative 
coupling as compared to the retarded regime (compare Figs. 7b 
and 4b). The non radiative coupling is dominated by electron scat- 
tering. For comparison purpose, we also represent the SPP contri- 
bution to the non radiative rate, even if the definition of the 
plasmonic mode becomes questionable in this regime (cf. Fig. 2). 
Due to strong losses, the 2D-LDOS no more follows a lorentzian 
profile as in Fig. 3b). Therefore, we estimate the SPP contribution 
using a lossless guide description since we demonstrated earlier that 
it correctly describes the coupling rate to the guided mode 26 . This is 
however an extrapolation from the propagating regime situation. It is 
nonetheless helpful for a qualitative discussion. 

The coupling rate into the SPP decreases very fast as a function of 
the distance (Fig. 7b). This abrupt decay reveals the strong mode 
confinement at the metal surface. We deduce from the slope the 
penetration of the mode in air to d SPP = 9 nm. This is confirmed 
by the calculated mode effective volume or area represented in 
Fig. 5b. In this quasi-static regime, the SPP mode volume reaches 
extremely low value V e ff ~ 10" 4 (/l/n 1 ) 3 similar to the localized SPP 
configuration 24 ' 38 . Nevertheless, strong losses limits the maximum 
achievable ^-factor to /? ~ 30% at d = 15 nm (not shown). 

Finally, Fig. 8 presents the cooperative decay rate for two quantum 
emitters near the metallic nanowire at the emission wavelength X = 
525 nm. Due to strong losses, SPP coupling does not exist anymore. 
We rather observe a low collective decay rate, governed by radiative 
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contribution since only direct free-space coupling occurs. Therefore 
the pseudo -periodic behavior of the cooperative decay rate follows 
the free- space wavelength. It is moreover strongly damped due to 
high losses in this regime. 

Discussion 

AID plasmonic waveguide efficiently promotes light-matter inter- 
action at the nanoscale. Although the SPP presents a modest quality 
factor (of the order of 10 to 100), they possess extremely low modal 
area (or volume) that are not diffraction limited and are thus very 
promising for enhancing single atom/single photon process. In the 
retarded regime, we demonstrated modal volume typically ten times 
below the diffraction limit and /? factor up to 90% but that could be 
still increased if crystalline nanostructures are involved. In this work, 
we also clarified the effect of Joule losses and proposed a simple 
physical understanding of the coupling mechanism in surface plas- 
mon coupled emission. We expect that this would be also helpful in 
defining strong coupling regime conditions 39 . 

Finally, we discussed the cooperative behavior of two coupled 
quantum emitters and demonstrated an efficient collective excitation 
of guided SPPs. It has been recently proposed to entangle two 
remotes qubits 1819 or to realize nano-optical logical gates 20 by con- 
trolling sub/superradiance. This regime is also of strong interest in 
the context of superradiant lasers (also called bad cavity lasers) pre- 
senting high spectral purity 40 ' 41 . Due to the poor quality factor of SPP 
resonances, we envision a fruitful analogy for designing plasmon 
nanolasers 42 44 . Indeed, the collective relaxation of several identical 
emitters in a poor cavity can lead to a coherent emission with extre- 
mely narrow spectral width 40 ' 41 . Practically, the limit linewidth of a 
single mode laser writes 45 



A kV N *P 

Av= - 

4nP 0Ut 



r 0 k 

To + k/2 



(4) 



where k is the cold-cavity loss rate, P out the laser output power, and 
N sp = N e /(N e — N g ) the degree of population inversion. In a good 
cavity regime (k<^2T), the laser spectral width is governed by cavity 

loss rate. It follows the Schawlow-Townes limit Av=—-^—. 

4tt P out 

However, as far as SPP amplication is concerned, one as to consid- 
ered the bad cavity regime (k^>2T): Av= —^JV — o jj^ s defines 

4tt P out 

the superradiant laser for which the electromagnetic coherence ori- 
ginates from the collective behaviour of the gain medium instead 
of the cavity feedback. Here, the cavity mode is only necessary to 



collectively couple the atoms forming the gain medium. The laser 
linewidth is then limited by the gain medium and could be largely 
below the Schawlow-Townes limit 40 ' 41 . Note also that a partial popu- 
lation inversion N sp degrades the linewidth lower limit. This has to be 
taken into account when describing SPP nanolasers since non radi- 
ative energy transfer to the metal limits the population inversion 46 . 
Since the SPP guided mode presents a rather low quality factor but 
mediates efficiently cooperative relaxation, we expect a fruitful ana- 
logy with superradiant lasers to achieve a coherent SPP nanosource. 

Last, we investigated the quasi-static regime and demonstrated 
SPP confinement down to 10~ 4 (/l/n 1 ) 3 but very low surface plasmon 
coupled emission due to critically important losses. Although limited 
to plasmonics nanowires in this work, the rich physics underlying the 
coupling between single photon source(s) and lossy materials is very 
general and remains valid for arbitrary geometries as well as meta- 
materials. Additionally, the extreme mode confinement observed 
at the band-edge of the dispersion relation is closely related to loca- 
lized surface plasmon and create a bridge with cavityless quantum 
electrodynamics. 

Methods 

Dispersion curve. The dispersion curve of the metal nanowire is achieved by solving 
the dispersion relation of the TM 0 i and HE n modes with Davidenko's algorithm in 
the complex plane 47 . This leads to complex propagation constant solutions 
k S pp(co) = k S pp + ik" S pp. The mode propagation length obeys L S pp = 1 /2k" S pp. 
Davidenko's algorithm is an efficient method to find complex root, including leaky 
modes (mode effective index below the surrounding optical index). In the present 
case, we observe that the HE n mode is not cut-off but its effective index approaches 
that of the surrounding optical index for small nanowire in agreement with the 2D- 
LDOS calculation (see Fig. 3) and also with the asymptotic expansion derived in ref. 8 
(see also supplementary information). 

Green's dyad. The 2D-Green's dyad is numerically computed discretizing the 
nanowire cross-section with rectangular meshes as detailed in ref. 48. Meshes sizes 
varies from 0.5 nm to 10 nm for larger radii. Importantly, the Green's dyad variation 
inside the mesh is exactly included, ensuring accurate evaluation. Arbitrary 
waveguide profile eventually supported on a wafer, are easily treated with this 
method 26 . Symmetrical properties of the Green's tensor are used to limit the 
numerical evaluation of G 2D on the k z > 0 range. Indeed, G 2D (r 1// , r 2// , — k z ) = 
[G 2D (r 2// , ry/, k z )] T (see supplementary information). The behaviour of G 2D near SPP 
resonance is also discussed in the supplementary information (Fig. S6). 

Decay rates. The dipolar total decay rate expresses T tot (x,y) /niT 0 = — 3/nik 0 
f ^Im [u. G 2D (r 0 //,r 0 //,/c z ) .u] dk z where r 0 // = (x, y) refers to the emitter position 
in the transverse plane and p = pu is the dipolar transition moment. The collective 
decay rate of two coupled dipoles (p x = piUi,p 2 = p 2 u 2 ) writes r^/wifo = —3/nik 0 

J ^Im [e lkzAz Ui-G 2D (r 1 / / ,r 2 / / ,k z ) ■ u 2 ] dk z with r v/ and r 2// the dipole positions in 
the transverse plane and Az = z x — z 2 their separation distance. Radiative rates T rad 
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and T™2 are obtained limiting the integration on the radiative waves (\k z \ < «i/c 0 ) an d 
cancelling the metal losses [Jm(e 2 ) =0]. SPP contributions T S pp and are 
achieved keeping the SPP resonance only into the integration. Thanks to the 
Lorentzian shape of the resonance, analytical expressions are available (see Eq. 2 and 
supplementary information). 

Relaxation channels in lossless or quasi-static limits. The Purcell factor near a non- 
absorbing waveguide is given by 21 ' 49 

T pU 3nce 0 \E u (R + d)\ 2 
T 0 ^ Aoo (ExH*).z.dA' 



(5) 



where (E, H) is the modal electromagnetic field, * refers to the conjugated form and u 
is the emitter orientation. For a circular cross-section, it is analytical (see 
supplementary information). 

In the quasi- static limit, radiative and non radiative rates expressions are given in 
ref. 29. 
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